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The role of active region coronal magnetic field in determining 
coronal mass ejection propagation direction 


Rui Wang!, Ying D. Liu!, Xinghua Dai??, Zhongwei Yang!, Chong Huang? and Huidong 
Hu! 


ABSTRACT 


We study the role of the coronal magnetic field configuration of an active re- 
gion in determining the propagation direction of a coronal mass ejection (CME). 
The CME occurred in the active region 11944 (S09W01) near the disk center on 
2014 January 7 and was associated with an X1.2 flare. A new CME reconstruc- 
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tion procedure based on a polarimetric technique is adopted, which shows that 
the CME changed its propagation direction by around 28? in latitude within 2.5 
Ro and 43? in longitude within 6.5 Ro with respect to the CME source region. 
This significant non-radial motion is consistent with the finding of Möstl et al. 
(2013). We use nonlinear force-free field (NLFFF) and potential field source sur- 
face (PFSS) extrapolation methods to determine the configurations of the coronal 
magnetic field. We also calculate the magnetic energy density distributions at 
different heights based on the extrapolations. Our results show that the active 
region coronal magnetic field has a strong influence on the CME propagation 
direction. This is consistent with the “channelling” by the active region coronal 
magnetic field itself, rather than deflection by nearby structures. These results 
indicate that the active region coronal magnetic field configuration has to be 
taken into account in order to determine CME propagation direction correctly. 


Subject headings: Sun: activity — Sun: coronal mass ejections — Sun: magnetic 
fields 
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1. INTRODUCTION 


Solar storms, known as coronal mass ejections (CMEs), are the most violent phenomenon 
in the solar atmosphere and a driver of major space weather effects. A reliable CME forecast 
has important consequences for life and technology on the Earth and in space. Accurate 
predictions of CME propagation direction are a major requirement for the forecast. However, 
CMEs do not always propagate in a constant direction and sometimes their propagation 
direction can be non-radial. CME propagation direction can deviate from the radial direction 
of the CME source region. Magnetic forces usually play an important role in deflecting 
CMEs within several solar radii. MacQueen et al. indicated_the influence of the 
background coronal magnetic field on CME deflection. Shen et al. and Gui et al. 
used gradients in the magnetic energy density of the solar corona to explain the 
observed deflection. Kay et al. developed a model called ForeCAT, which 
considers the influence of magnetic forces on CME propagation. Panasenco et al. (2013 
demonstrated that asymmetric large scale magnetic forces acting on CMEs can make CMEs 
propagate in the direction of the least resistance. The deflections can occur in both latitude 
and longitude '2010alb 5010. 013). Isavnin et al. 
reconstructed 15 CMEs between the minimum of Solar Cycle 23 and the maximum 
of Solar Cycle 24. They found that latitudinal deflections up 35? are observed while the 
maximum longitudinal deflection is only 5?.4. Without out-of-ecliptic observations, more 
latitudinal deflection events were reported than longitudinal deflections 
198. tpa et all koos sane et all 2018 [Panasenco otal] 2012). Byme et al. (oio 
show a latitudinal deflection of —30? below 7 R for the 2008 December 12 CME. With 
out-of-ecliptic observations provided by the Solar Orbiter satellite that will be launched into 
space in the near future, longitudinal deflection events will gain more attention. 


According to previous studies, coronal holes where the fast solar wind flow comes from 
are thought to be one of the important factors that can deflect CMEs lCremades et_alll2008). 
The coronal hole acts as a magnetic wall that constrains CME propagation. Kilpua et al. 
suggested that CMEs cannot penetrate the magnetic wall and are guided by polar 
coronal hole fields to the equator. Gopalswamy et al. (2009 considered the influence of 
the magnetic field strength of coronal holes and defined a coronal hole influence parameter 
(CHIP), whose direction has a good agreement with the direction of CME deflection. On the 
other hand, the CME moves away from nearby coronal holes and tends to deflect towards 
and propagate along a nearby helmet streamer or streamer belt (SB), for there are X-type 
magnetic structures (X-lines) where the magnetic stresses are minimum peni enn 
TORE STI TETI. Xie et al. (2009) found that slow CMEs (speed less 
than 400 km s^!) tend to propagate towards and then along the SB due to the deflections 
by strong polar magnetic fields of corona holes. Kay et al. (2013 suggested that magnetic 
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background that contains strong magnetic gradients deflects CMEs to reach the SB; however 
for a weaker background, it just moves CMEs towards the SB but may not be capable of 
fully deflecting the CME along the SB. Panasenco et al. and Kay et al. both 
think about the CME deflection away from coronal holes and towards a streamer. Panasenco 
et al. (2013 also looked at the rolling motion of filaments, which is similar to the deflection 
of CMEs. 


As mentioned above, the driving forces originate from the magnetic field of the coronal 
holes outside the active region (AR), i.e., the CME source region. However, the role of 
the magnetic forces coming from the AR itself is not well understood. AR NOAA 11944 
(S09W01) produced a fast halo CME with an X1.2 flare (S12W08) on 2014 January 7. This 
is an important event from a space weather perspective as forecasters expected a significant 
impact (a G3 class geomagnetic storm or higher) on the Earth, yet it did not occur. Ac- 
cording to the study of Möstl et al. (2015), the AR was near the disk center, but the CME 
propagated towards Mars, which was 51° away from the Earth in longitude. The deviation 


angle of the CME non-radial propagation in longitude was much larger than we thought. 
Most cases that have been reported have coronal longitudinal deflections <~ 20° 
boni Kav et allor Isovnin et all aor slthough, there may be coun 
terexamples in observations and simulations. Most] et al. (2015 showed that the distance of 


the coronal holes to the AR is large, which gives a CHIP value (defined by Mohamed et al. 
2012) smaller than the one necessary to deflect a CME almost completely away from the 


Earth. They used a Graduate Cylindrical Shell (GCS) model (Thernisien et al! and 
a triangulation technique based on the STEREO data to estimate the 
CME propagation direction. Their results reveal that the change in CME propagation direc- 
tion occurred close to the Sun and was likely to be related to the AR magnetic field rather 
than the coronal hole. They presented some global evidence to prove that the coronal mag- 
netic field “channelled” the CME rather than deflected the CME. We follow their definition 
of the coronal “channelling”, which is non-radial starting already with its inception on the 
Sun due to the AR magnetic fields. The “deflection” implies a change in direction, caused 
by the nearby structures, e.g., coronal holes and streamers. Here we present local evidence 
around the AR related to the coronal “channelling” as well as new CME reconstruction 
results. The topology of the magnetic field at different altitudes, especially near the CME 
source region, and the influence of magnetic stresses of the strong AR magnetic fields are 
both worth further studying. 


In this paper, we use a polarimetric technique to reconstruct the CME and compare the 
results with those of Mostl et al. (2015). We also investigate the configuration of coronal 
magnetic fields near the source region with a nonlinear force-free field (NLFFF) method 
and the magnetic field topology on a larger scale using the potential field source surface 


e. 


(PFSS) technique. Meanwhile, the changes of magnetic field energy density distribution with 
heights will be examined, which can demonstrate how magnetic forces play a significant role 
in channelling the CME. 


2. OBSERVATIONS AND ANALYSIS 


AR 11944 (S09W01) produced a very fast halo CME (linear speed of ~1830 km s^!) 
and an X1.2 flare (S12W08) on 2014 January 7. The flare started at 18:04 UT, peaked at 
18:32 UT and ended at 18:58 UT. There exists another AR 11943 (S11W19) southwest of AR. 
11944. The CME source region is located around the border of the two ARs and belongs to 
AR 11944. At the time of the event, the two STEREO spacecraft were 151° ahead (A) and 
153? behind (B) in heliospheric longitude with respect to the Earth and at distances of 0.96 
AU (A) and 1.08 AU (B) from the Sun, respectively. Figure 1 shows a sided view of the CME 
from COR1 of STEREO-A and a backsided view from COR1 of STEREO-B. Note that the 
AR 11944 is located near the disk center as observed by the Atmospheric Imaging Assembl 

AIA; on board the Solar Dynamics Observatory (SDO; 
"mi We can see that the CME was propagating to the southeast direction in STEREO- 
B. The non-radial propagation of the CME in longitude is more clearly shown in COR2 of 
STEREO-B where the CME was propagating to further distance. From the viewpoint of the 
Earth, the CME would be propagating to the southwest. Most] et al. demonstrated 
that the CME propagated towards Mars, which was at a heliospheric longitude of 51? west 
of the Earth. They showed that the CME was deflected by 37+10° west of the source region 
or 45+10° away from the Sun-Earth line in heliospheric longitude (Masti et_all{2015). They 
further pointed out that the final direction of the CME was attained very close to the Sun, 
around 2.0 R from the center of the Sun. 


Figure 2a shows the EUV structures of the CME source region before the eruption 
occurred. Helioseismic and Magnetic Imager (HMI; on board the SDO 
provides us high time and spatial resolution magnetograms with 0".5 per pixel, as shown by 
the contour map overlaid on the EUV images. The positive polarity sunspot in the left side of 
Figure 2a has a very strong photospheric magnetic field, whose vertical component exceeds 3 
kG. Livingston et al. (2006 found that the percentage of 12,804 sunspot groups between 1917 
and 1974 with magnetic fields 23 kG is only 4.696. Some open EUV structures connected 
to the positive polarity sunspot are outlined by the red dotted lines and two lower EUV 
loops are outlined by the blue dotted lines. Figure 2b-2d show successive running difference 
images during the CME eruption, which better exhibit the evolution of the features in the 
kernel of the CME source region. The global solar disk observations of the erupting CME 
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have been examined by Möstl et al. (2015). The flare, coronal dimmings, post-eruption 
arcades, and a global coronal wave are thought to be the classical disk signatures of an 
erupting CME. The location of the front of the global coronal wave that is thought to be 
driven by the lateral expansion of the CME is visualized in their work. The global coronal 
wave at the onset phase can also be detected through the evolution of the features in Figure 
2b-2d. The lower EUV loops in Figure 2a were blown away by the coronal wave propagating 
towards southwest. Between the strongest magnetic fields (red contours) at AR 11944 and 
the negative photospheric magnetic fields (cyan contours) at AR. 11943, a possible magnetic 
flux rope (MFR) rose up then erupted with the wave, and was followed by a flare (see Figure 
2d). In the west of the positive photospheric magnetic field (white), the coronal wave front 
is outlined by the red dashed line (see Figure 2c). The disk wave front likely corresponds to 
the CME front. According to the propagating direction of the wave, we can estimate that 
the CME was propagating away from the positive polarity sunspot with strong magnetic 
fields at the initial stage of the eruption. This is consistent with the coronagraph image in 
Figure 1b. Based on the evolution of the features above, we can roughly confirm the location 
of the CME source region as the same position of the flare. 


2.1. POLARIMETRIC RECONSTRUCTION OF THE CME 


Here we use a polarimetric technique to reconstruct the CME. This technique was devel- 
oped by Moran et al. to reconstruct the 3D density distribution from the combination 
of intensity and polarized brightness images of LASCO. Based on the Thomson scattering 
theory, the polarization ratio determines the line of sight (LOS) averaged distances from 
the plane of sky (POS), which corresponds to the polarized image of CMEs as a reference 
surface for 3D polarimetric reconstruction. According to the reconstructed 3D location, each 
reconstructed CME pixel in the view of the equatorial plane can be located. The mass at 
each pixel depends on the position of each reconstructed pixel and the total brightness of the 
corresponding pixel. The mass at each pixel can be calculated using Equation (9) of Dai et 


(2015). 


; 2015). The advantage of the polarimetric reconstruction method is that, by using 


al. The technique has been further improved for mass calculation of CMEs by Dai et 


al. 
polarized observations from only a single viewpoint, it can reconstruct the mass distributions 
of CMEs. However, there exists ambiguities in the polarimetric analysis of CME observations 


(Dai et at]2014). Not all data can be utilized for reconstruction. For example, coronagraph 
images from STEREO-A show an implicit ambiguity (Dai et a1] 014), and thus cannot be 
used. Figure 3a shows that the CME is split into two parts by the POS of STEREO-A (red 
dashed line). In the polarimetric technique, we cannot distinguish the sign of a in Equation 
(15) of Dai et al. , which is the angle between the POS and a radial line connecting 


Table 1: CME reconstruction results at different times. The longitude is counted from the 


Sun-Earth line. 
Mass center Time (UT) Mass of CME (g) Heliocentric Longitude (°) Latitude (°) 


\Instrument distance (Ro) 
P1 18:25\B-COR1 1.08e+15 2.07 22 -37 
P2 18:30\B-COR1 2.68e4-15 2.24 22 -39 
P3 18:35\B-COR1 4.06e+15 2.34 25 -40 
P4 18:40NB-COR1 5.11e4-15 2.29 28 -39 
P5 18:45\B-COR1 5.87e+15 2.33 29 -38 
P6 20:00\B-COR2 1.10e+16 6.54 51 -40 


the center of the Sun and the average 3D location of electrons along the LOS of the CME 
pixel. Because of cos(-a)=cos(a), a € [0°, 90°), the sign of a is ambiguous. The ambiguous 
locations of the electrons along the LOS of the CME pixel are denoted as the front and back 
of the POS corresponding to œ and -a. Thus we have to assume that the entire CME is 
either located before or behind the red dashed line. Therefore, although COR1 gives a sided 
view of the CME as shown in Figure 1a, it cannot be used for the polarimetric reconstruction 
due to this ambiguity. Fortunately, the viewpoint from STEREO-B has no such ambiguity 
because the entire CME is behind the POS of STEREO-B (see Figure 1), and thus it is 
reliable and appropriate for reconstruction. 


The mass centers of the CME reconstructed from COR1 of STEREO-B at different 
times are also shown in Figure 3a. The green solid line represents the Sun-Earth line. The 
red and blue lines represent the LOS of STEREO-A and B, respectively. The biggest solid 
black sphere represents the mass center of the CME at the earliest time, i.e., 18:25 UT. The 
yellow lines are the radial lines connecting the mass centers with the center of the Sun. Table 
1 gives the results of our reconstruction using data from COR1 of STEREO-B. Within 2.5 
Ro, COR1 of STEREO-B provides us successive polarized data every 5 minutes. We use 
COR1 data because it is more accurate for its smaller field of view (FOV) (as shown in Figure 
1) than COR2. After 18:45 UT, the CME exceeded the FOV of COR1 of STEREO-B, so we 
have to consider COR2 data. However, COR2 of STEREO-B cannot provide us polarized 
data until 20:00 UT (as shown in Figure 1c). After 20:00 UT, the CME exceeded the FOV 
of COR2. Therefore, for COR2 of STEREO-B only the polarized data around 20:00 UT 
can be used. From 18:25 UT to 18:45 UT, the mass centers of the CME are at the heights 
between 2.0-2.5 R relative to the center of the Sun (see Figures 3a). The mass of the CME 
keeps increasing with distance as expected. The heliospheric longitude of the CME mass 
center is slightly increasing but the latitude keeps around 40? (south). The CME deflection 
in heliospheric latitude can be seen from the coronagraph images of STEREO-B in Figure 
1. Our results suggest that the CME was already deflected to a large longitude and latitude 
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within 2.0 Rg. Compared with the approximate CME initial place (S12W08), the CME was 
deflected in heliospheric longitude by more than 20° below 2.5 Reo. Moreover, around 6.5 
Ro the longitudinal deflection reached around 40° with respect to the source region (also 
see Figures 3a and 3b). As shown by previous studies, most coronal longitudinal deflection 
events generally have a deflection <~20°. The longitudinal deflection in the present case is 
unusually large in CME deflection observations. 


Figure 3c and 3d show that the full CME white light structure in Figure 1c is well 
reconstructed by our polarimetric technique. Figure 3c shows the distribution of the recon- 
structed CME mass from COR2 of STEREO-B view. Figure 3d shows the reconstructed 
CME pixels in different heliocentric distance denoted by the color bar. The black solid poly- 
gon shows the boundary used for determining the CME mass center. When we determine 
the boundary above, it will introduce human error, but compared with the systematic un- 
certainty it is negligible. It is difficult to determine the systematic uncertainties inherent in 
our reconstruction method. However, our experience using the method suggests an error of 
about 10° for both the longitude and latitude. The results of Mostl et al. (2015 show that 
the CME is attained to W32+10° between 2.1 and 18.5 R by means of the GCS method, 
and attained to W45+10° between 20 and 30 R by means of the triangulation method. Our 
results are consistent with theirs below 2.5 Ro, but the CME has already been channelled 
into a direction around W50° at 6.54 Ro (see Table 1). The reason is that the polarimetric 
method uses the mass center to show the position of a CME rather than the CME apex used 
by the GCS method. From Figure 3c we can see that a fraction of the mass is concentrated 
at the northeast of the Sun, which may result in the deviation from the results of Most] et al. 
(2013). For the reconstruction from COR1, it does not encounter this problem. On the other 
hand, Table 1 shows that the CME maintained around $40? in latitude beyond ~2.0 R but 
kept moving westward in longitude. This may be related to the large scale configurations 
of the magnetic field. The streamer belt looks like a ^U-shape" structure surrounding the 
source region and opens to the west. The CME moved to the south until reaching the top 
of the streamer belt (the X-line), then moved along the X-line. Therefore, the latitude kept 
constant beyond ~2.0 Ro. 


Our reconstruction results are consistent with those obtained by Most] et al. (2015), who 
used a forward modeling technique and a triangulation method. According to the synthesis 
of heliospheric observations of Most] et al. (2015 , the CME arrived at Mars which was at 
a heliospheric longitude of 51° west of the Earth. From our CME reconstruction, we also 
see that the CME has already been channelled to a large angle at a relative low altitude. 
We think that the strong coronal magnetic field from the positive polarity sunspot is likely 
to play a role in determining the CME propagation direction. In the following section, 
we will investigate the configurations of the coronal magnetic field using two extrapolation 


techniques. 


2.2. CORONAL MAGNETIC FIELD RECONSTRUCTION 


We use different extrapolation models to reconstruct the coronal magnetic field within 
and surrounding the ARs. The simplest way to model the coronal magnetic field is to 
assume that the field is potential, i.e., no electric current. An AR, however, usually contains 
a large amount of electrical currents. The nonlinear force-free field (NLFFF) method is 
more appropriate in determining the coronal magnetic field configuration for an AR, while 
the PFSS technique may be more appropriate for a global field reconstruction. 


Within ~0.16 Ra from the photosphere, we adopt the NLFFF method as proposed b 
and extended by and 
. We use the NLFFF optimization code, which has been optimized for the application 
to SDO/HMI vector magnetograms [Wiegelmann et at] 012), to extrapolate the coronal 
magnetic field from the observed vector magnetograms in a Cartesian domain. A prepro- 


cessing procedure (Wiegelmann et al14|2006) is employed to remove most of the net force and 


torque from the data, so the boundary can be more consistent with the force-free assump- 
tion. This extrapolation technique has been successfully applied in our previous work on 
two successive solar eruptions nee: abou. Here, we apply this NLFFF method to an 
area of 700" x 360". The spatial resolution is 1".0 per pixel. We adopt the Lambert Cylindri- 
cal Equal Area projected and remapped vector magnetic field data I ns oe 
[Calabretta & Greise 2003 from HMI 


Figure 4 shows the NLFFF extrapolation results covering the main source region from 
different viewpoints. In order to distinguish different magnetic flux bundles, we use different 
colors to represent them (see Figure 4a). Compared with the EUV images of Figure 2b, the 
open magnetic flux bundles in black are consistent with the open EUV structures outlined 
in Figure 2b. The flux bundles in white are the core structure, which is shown as a fan-like 
structure inclining to the west. The flux bundles in blue are consistent with the post flare 
arch structures. The source region location is shown as the yellow star. As can be seen, 
the source region is located between the positive polarity sunspot and the open fan lines. 
The CME was probably affected by the strong magnetic pressure from the positive polarity 
sunspot and inclined to propagate along the open fan structures. However, we do not get 
the possible flux rope structure appearing in Figure 2a. A possible explanation is that this 
structure is very dynamic and therefore cannot be reconstructed with a model based on a 
static snapshot. In fact, we find some twisted field lines entangled with the lines originating 
from the sunspot with strong magnetic fields, which may evolve into this structure eventually. 


—9Q- 


From 1.16 to 2.5 Ra, the potential field source surface (PFSS) method 
1969: is applied to reconstruct the coronal magnetic topology in 
a larger scale. Figure 5a gives the large scale structures of the magnetic field surrounding 
the ARs. The open magnetic field (cyan) is surrounded by the SB (purple). Figure 5b shows 
the open magnetic structures and some closed magnetic structures as part of the SB below 
2.5 Ro. The open magnetic structures are similar to the NLFFF results at higher altitudes. 
One end of the open magnetic field lines roots in the strong positive magnetic field region, 
and the other end roots near the fan structure. The closed field lines are above the ARs as 
part of the SB and connect the positive polarity sunspot with strong photospheric magnetic 
fields to the negative polarity sunspot in the east of the AR. 11944 (see Figure 4a). From 
Figure 5b we can see that the open magnetic structures are not along the radial direction 
but incline to the west. Figure 5c shows that the open field lines stretch out towards the 
north and south poles, respectively. The bundles of the open field lines stretching out to 
the south just pass through our reconstructed CME mass centers, which implies that the 
trajectory of the CME was likely along the open magnetic fields originating from the ARs. 


The Lorentz force comprise a pressure force and a tension force. The tension force is 
the means to trap coronal plasma against the expansion of CMEs. As Figure 4b and c have 
shown, the closed magnetic field lines containing the tension force cover the source region 
and prevent the CME from expanding to higher places. The region of the open field lines, 
however, creates much less resistance, so the CME is more likely to pass through there. 


Previous work (Xie et al 009: Zuccarello et al 2012; [Kay et al] 013) indicated that 


the CME moves away from nearby coronal holes and tends to deflect towards and propagate 
along a nearby SB. For the distribution of the SB of this event, the CME propagated towards 
the SB in latitude but moved away from it in longitude (as shown in Figure 5a). According 
to Xie et al. (2009), there is no obvious correlation between the deflection of fast CMEs (V > 
400 km s^!) and the SB. The remarkable difference of this CME deflection event from others 
is that there exists a sunspot with unusually strong magnetic fields nearby the erupting CME. 
This would produce a strong magnetic barrier and force the CME to propagate away from it. 
In order to evaluate the influence of the strong magnetic field on the CME propagation, we 
investigate the magnetic energy density distribution with altitudes in the following section. 


2.3. MAGNETIC FIELD ENERGY DENSITY DISTRIBUTION 


Previous studies showed that a nonuniform distribution of the background magnetic 
field energy density may cause the deflections of CMEs at an early stage and make CMEs 
tend to propagate to the region with lower magnetic energy density. Panasenco et al. (2013 
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used magnetic stresses to represent this. The model of Kay et al. for forecasting a 
CME’s altered trajectory adopted magnetic pressure gradient. For the sake of simplicity, 
we present the magnetic energy density at some typical heights to study the role of the AR 
magnetic field in forcing the CME to propagate in a non-radial direction. The magnetic 
energy density Dme is expressed as follows, 


B? 


Das = a 
20 


Uo = 41 x 10 7W/A.m 1 (1) 


where B is in spherical coordinates, i.e., B? = By? + Bs + B. 


Figure 6 shows the magnetic energy density at 1.06 and 1.16 Ro, respectively. At such 
heights, the magnetic energy density comes from the results of the NLFFF extrapolation. 
The CME source region (812W08) is located at the southwest of the positive polarity sunspot 
with strong magnetic fields (as shown in Figure 2). The distribution of the magnetic energy 
density shows a very strong gradient towards the southwest. This is produced by the presence 
of the strong sunspot. With the height increasing, the two-peak energy density distribution 
above the AR 11944 in Figure 6a gets into one peak in Figure 6b. The projections of 
the footpoints of the magnetic field lines traced from the CME mass centers (whose colors 
correspond to those in Figure 5) are located at the AR 11943. 


Figure 7 gives the magnetic energy density distributions at higher altitudes, based on 
the PFSS extrapolation. Combining Figures 6 and 7, we find that the relative positions 
between the peaks of the density distributions (white contours) and the source region (the 
black cross) changed with the height. The magnetic field from the strong field region acts 
as a magnetic wall that constrains the CME propagation. Both Figures 6 and 7a show 
that the gradient descent direction points to the southwest, which is consistent with the 
CME propagation direction. Figure 7b shows that there appears another density peak in 
the northeast at 2.0 Ro. The gradient descent direction of the energy density distribution 
near the CME mass centers gradually changed from the south to the west. This gradient 
is consistent with the change of the projections of our reconstructed CME mass centers 
(from the colored crosses to the white one). Therefore, we think that the change in the 
CME propagation direction beyond ~2.5 Reo was likely caused by the presence of second 
density peak, while the change in the CME propagation direction within ~2.5 R probably 
resulted from the positive polarity sunspot with strong magnetic fields. The location of these 
open fields corresponds to the region of the lowest magnetic energy. The direction from the 
first density peak to the density valley is consistent with the propagation direction of the 
CME (from the black cross to the colored crosses; see Figure 7), which implies that the 
CME propagates to the less resistance region. Sum up, the change of the relative positions 
between the CME source region and the projections of the reconstructed CME mass centers 
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is consistent with the change of distributions of the magnetic energy density with heights. 
This consistency supports the idea that the non-uniform distribution of the coronal magnetic 
field around the ARs plays a key role in channelling the CME. Note that we do not take into 
account any change in the magnetic field topology due to the passage of the CME. However, 
this method may still identify the environment that has influenced the CMEs channelling. 


3. SUMMARY AND CONCLUSIONS 


In this work, we concentrate on the role of the active region coronal magnetic field in 
channelling a CME. A new CME reconstruction procedure based on a polarimetric technique 
is adopted. Our polarimetric reconstruction results imply that the CME was deflected close 
to the Sun and attained to its final direction within several solar radii. Our results also show 
that the CME was deflected in both latitude and longitude. The CME is demonstrated to be 
channelled into a direction around 28° (in latitude) away from the source region within 2.5 
Ra. In longitude, it was channelled into a direction around 20° with respect to the source 
region within 2.5 Ro. This longitudinal deviation become about 40° from 2.5 to 6.5 Ro. 


We use the NLFFF and PFSS methods to extrapolate the coronal magnetic field around 
the active regions at different heights. The NLFFF extrapolation results show a fan-like open 
magnetic field structure in the west of the source region. We obtained similar open field 
structures in higher corona from the PFSS model. We think that these open field structures 
probably provided a favorable channel for the CME propagation. In order to evaluate the 
influence of the magnetic pressure on the CME propagation, we examine the magnetic en- 
ergy density distributions at different heights based on the NLFFF and PFSS extrapolation 
results. We compare the corresponding positions of the magnetic energy density peaks of 
each layer and the CME, as well as the relations between the gradient descent direction of 
the energy distribution and the changes of the CME mass centers. These comparisons indi- 
cate that the positive polarity sunspot on the northeast of the source region, which shows 
photospheric magnetic fields larger than 3 kG, plays a significant role in determining the 
CME propagation. There are two main factors that may have caused the CME non-radial 
propagation: (1) the magnetic pressure from the strong positive polarity sunspot forced the 
CME to the southwest; and (2) the configurations of the open and closed magnetic fields 
around the source region made the CME propagate to the least-resistance direction. 


To conclude, this study shows the important role of the active region coronal magnetic 
field itself in determining the CME propagation direction. The active region coronal magnetic 
fields cause a significant change of CME propagation direction in both latitude and longitude. 
While our results are consistent with the “channelled” CME motion defined by Most] et al. 
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2015), they are obtained by a new approach and give a more complete story that is necessary 
to understand how the active region coronal magnetic field has channelled the CME at low 
altitudes. These results will help determine the propagation direction of CMEs correctly, 
which is important for space weather forecasting. 


The work was supported by the Specialized Research Fund for State Key Laboratories 
of China, NSFC under grant 41374173, and the Recruitment program of Global Experts of 
China. 
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Fig. 1.— Running difference images of the CME from STEREO-A and B. (a) A sided view 
of the CME from CORI of STEREO-A. (b) A backsided view of the CME from COR1 of 
STEREO-B. (c) COR2 view of the CME from STEREO-B. The CME was propagating in 
the southeast direction for STEREO-B. The white circles mark the location and size of the 
solar disk. 


—16- 


an erupting MFR? 


0 

-5 
mS m 
o0 o0 
[^] oO 
= S 
o —10 v 
"d "B 
d d 
~ ~ 
= em 
S -15 S 

—20 


5 10 15 20 25 0 10 20 30 
Longitude (deg) Longitude (deg) 


Latitude (deg) 
Latitude (deg) 


0 10 20 30 0 10 20 30 
Longitude (deg) Longitude (deg) 


Fig. 2.— The EUV structures and eruption of the CME source region from AIA 171 A. (a) 
EUV structures from AIA 171 A. The red and blue dotted lines mark open EUV structures 
and short EUV loops, respectively. (b)-(d) Successive running difference images of the 
eruption. The red dashed line marks the wave front of the CME. Contours of -1500, - 
500, 500, 1500, and 2500 G of the photospheric magnetic field along the line of sight from 
HMI magnetograms are overlaid on the EUV images and are marked in purple, cyan, black 


(a) (white (b)-(d)), blue and red, respectively. 
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Fig. 3.— Reconstruction of the 2014 January 7 CME from STEREO-B images. (a) Recon- 
structed CME mass distribution from COR1 of STEREO-B viewed from solar north. Green, 
red, and blue solid lines indicate the directions of the Earth, STEREO-A, and STEREO-B, 
respectively. The black solid and yellow lines show the directions of the mass centers (five 
black spheres) of the CME from 18:25 UT to 18:45 UT, respectively. The mass center at 
the earliest time is shown by the biggest black sphere. The dashed colored lines denote 
the corresponding plane of sky. The reconstruction results (color shading) at 18:25 UT are 
used. (b) Reconstructed CME mass distribution from COR2 of STEREO-B viewed from the 
Earth. The reconstruction results at 20:00 UT are used. The 171 Å image is mapped on the 
solar surface. The mass of each CME pixel is indicated by the color bar. (c) The distribution 
of the reconstructed CME mass from COR2 of STEREO-B view. (d) The distribution of 
the reconstructed CME pixels at different heliocentric distance (indicated by the color bar) 
from COR2 of STEREO-B view with error estimation as shown by the black solid polygon. 
The Sun is represented by the small red circle. 
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Fig. 4.— NLFFF reconstruction of coronal magnetic fields of the source region at 18:00 UT. 
(a) Different magnetic flux bundles are represented by different color field lines. The yellow 
star marks the location of the CME source region. Color bar gives the vertical component 
of the background magnetic fields of the photosphere. (b) and (c) are viewed along the 
directions shown by the blue and red arrows in (a), respectively. 
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Fig. 5.— PFSS reconstruction of the magnetic field above the source region. (a) Open 
magnetic fields (cyan) surrounded by the SB (purple) as viewed from the Earth. The observed 
magnetogram is mapped on the solar surface. The CME mass centers are represented by 
green, golden, red, cyan, and purple spheres, respectively. (b) Open magnetic field lines 
viewed from the south. Some closed magnetic field lines as part of the SB are also shown. 
(c) The same structures as shown in (b) but viewed from the west. 
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Fig. 6.— Distribution of the magnetic energy density at the heights of 1.06 (upper panel) 
and 1.16 Rọ (lower panel) based on the NLFFF reconstruction. The black cross shows 
the projected location of the CME source region. The colored crosses (P1-P5) show the 
projections of the footpoints of the magnetic field lines traced from the reconstructed CME 
mass centers. 
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Fig. 7.— Distribution of the magnetic energy density at the heights of 1.5 (upper panel) 
and 2.0 R (lower panel) based on the PFSS reconstruction. The black cross shows the 
projected location of the source region. The colored (P1-P5) and white (P6) crosses show 
the projections of the reconstructed CME mass centers. 


